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We investigate the problem of self-gravitating dust disk dynamics in a static state taking into 
account nonlinear effects. For this purpose Schrodinger-type equation including the mass conser- 
vation law is used for the whole description of hydrodynamic flows of self-gravitating dust. We 
have shown a purely hydrodynamic mechanism of ring formation in the radial direction by taking 
into account nonlinearity in the lowest order of expansion parameter, which determines an order of 
magnitude of flow. 



1 Introduction 

Self-gravitating systems are of great interest for investigation in astrophysics 
because of their widespread appearance P El E|- Such systems are difficult 
for analytical analysis when we consider a lot of factors, and therefore the 
standard way of analysis is in the terms of density perturbations. As a result 
the problem may be entirely linearized, and it facilitates analysis, but in this 
way one cannot entirely investigate some specific effects caused by nonlinear 
dynamics of such systems. In this paper the problem is also considered in 
terms of perturbation theory, but it is possible to take into account nonlinear 
effects of dynamics in the lowest order by means of elimination of secular in 
time expansion terms. Schr5dinger-type equation^ representing combined de- 
scription of Euler fluid dynamics together with the continuity equation is used 

^Of course we consider this equation only as auxiliary one without Plank's constant and any "quantum 
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for determination of potential hydrodynamic flow. Dynamics of rather outly- 
ing disk regions is discussed. These regions are situated far from the central 
massive object and the disk surface. Internal boundaries of these regions are 
predetermined by congruence condition of the self-gravitating dust potential 
and the compact body potential. The main aim of this paper is to consider a 
ring structure formation due to nonlinear hydrodynamic flow of self-gravitating 
dust. The internal region dynamics nearby the central object and the disk sur- 
face must be considered as the internal solution problem similarly to boundary 
layers problem. 



2 The usage of Schrodinger-type equation in Euler fluid 
dynamics 

Let us consider Schrodinger-type equation: 

i^t + a(^) A* - [/(x, i)^ = 0. (1) 

Here is a dimensionless complex wave function, i is an imaginary unit, a = 
a{t) is a certain dimensionless real function of time U (x, t) is a real potential- 
like energy function, A is three-dimensional Laplacian. 

Besides the ([T]) it is necessary to consider the complex conjugate equation: 

-z^* + a{t)A^' - t/(x, t)** = 0. 

Multiplying the previous equation by ^ and ([T]) by and subtracting the 
second expression from the first one, we get 

+ div (zal^pVlnG) = 0, (2) 



where G = "^*/"^. Dividing dH) and the conjugate equation by and ^* 
correspondingly, and summing these expressions gives 

d ^ A* A** 
i— InG + a—— + a— 2U = 0. 3 
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can be interpreted as differential conservation law for density l^'p of fluid 
moving with the velocity v = iaVlnO. 

For an arbitrary Euler flow v we have the following identity: 



vt + (v, V)v = ^V|v|^ + [v X rotv] + Vf. 



(4) 



Note that in our case the field v is potential: rotv = 0. For analysis of right- 
hand side of (jH) one can use (EI). As a result it follows the identity 



1^1 ,2 



+ 2aV[/ + -V. 



This result through identity dH) can be combined into the following Euler's 
equation for the flow v 



vt + (v, V)v = aSJ 



(5) 



where n = —ajoi = K{i) is a coefficient of a linear friction as it is called in fluid 
dynamics. Corresponding frictional force in the system —k\ can be regarded 
as a result of collisions between dust particles. In fact, this rather simplistic 
approach doesn't give an exact picture. Nevertheless, it allows to estimate the 
influence of a dissipation on a stationary dust distribution. Below we discuss 
the case = /^o = const. Dependence Av on time can be interpreted as particles 
interaction to be changed in time, e. g. due to increasing particles in size . 

Consider force per unit mass in the momentum equation in details. One can 
see from (EI) the force is potential. In real fluid dynamics, however, right hand 
side in Euler equation in the presence of the field of a potential force is the 
following 

F = --VP-V0, (6) 
P 

where p is density of fluid or gas , P is pressure, is potential of force per unit 
mass is considered below to be Newtonian force. We investigate dust objects in 
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this paper. And the state equation for dust is well known as p = 0. Therefore, 
force per unit mass (jHI) consists of dust self-gravitation and it may include 
gravitation force from a massive object, which is nearby the dust. Then we 
find the relation between Euler fluid dynamics and Schrodinger-type equation 

- 2a^^^ + 2aU = -cj). (7) 

3 The hydrodynamic equations for self-gravitating dust 

The problem considered in this paper can be formulated in the following way. 
We investigate dust objects having the equation of state P = 0. And dust is in 
self-gravitation in terms of Euler fluid dynamics. So, we come to the system of 
equations for self-gravitating objects dynamics, which consists of Schrddinger- 
type equation dH), the concordance equation and the Poisson one, which is: 

A0 = 47tG{po\^\^ + a6{z) + Mo^(r)), (8) 
where po - a characteristic density such that the function of density is 

The second z-direction ^-like source item in right hand side of the equation 
for potential describes a matter originally concentrated in thin disk with surface 
density a = a{x, t), where x, y are Cartesian coordinates at the disk surface. 

Now we shall make a nondimensionalization as following: 

f = r/Po, r = t/To, = 0/00- 
And make a notion : Qf(r) = aof{r), where 

T 

f{r) = exp{- J K{r)dr'} 
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is a dimensionless function of time. The system of equations then reads 

itr + £/(r)A* -WH = 0, (9) 

-2-'V^{r)^ + ^'f(-W = (10) 
A* = ^^^|*P, (11) 

00 

where e = Toao/^o " dimensionless characteristic parameter, estimating 

an order of magnitude of the dimensionless flow in the system: 

W = UTq is the non-dimensional collective potential. 

Below we omit ~, implying the equations to be written in a dimensionless 
form. 

In this paper we are interested in cases such that e << 1 is a small param- 
eter, i.e. the velocity of flow is small and the system in a gravitation field is 
about equilibrium . 

Researching of shows that for the equations to describe non-trivial 

situation we must use the conditions as follows: 

2 A^n^ d2 



Kq 00 

where /i is the first-order constant to e. From this it follows 0o = 
eRI/Tq, pq = £ fi / (AttGTq) , i.e. the dust density and the self gravitational 
potential are small and have the same order with respect to e. 



4 Approximate equations 

Let us seek the solutions as power series in e: 



00 



n—\ n=l n=l 
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Since small parameter e in (El) and f fTOl l is at the derivative of higher order, 
we can expect boundary layers in the system to appear. They could exist in 
the center of the field and at the surface of the disk. This boundary layers are 
connected with nonlinear mode not with viscosity. Outside of this boundary 
layers, i.e. far from the field center and the disk surface we use ordinary axial 
coordinate z. Nearby the disk surface we must use coordinate Z = z/e. Note, 
(5-like source in ( ITTl ) must be taken in account only in internal solution. 

For external region we have the system of equations by substituting the 
expansions in equations in two first orders : 

z*o,r = VFo^o; 2/(r)iyo = -^o; A$o = /i|*o|'; (12) 

-2/'(r)^J[^^ + 2/(r)[/i = $i, A$i = + *t*o). (13) 

l^ol 

Suppose the flow and the gravitational field in lowest order are stationary; then 
we have solution in this order: 

^0 = C'o(r) exp |-z J Wo(r, r)dr} , Wq = - ^^^O' 

where function $o is to be obtained from the Poisson equation: 

A$o = /i|Cop. (14) 

After some simple manipulations, we arrive at the following solutions at the 
first order : 



r 
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2/(r) ' ' ' \Co\ 
Here 
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Substituting these expressions in the equation for $i (fT3l)we get: 

A$i = /i(CoC* + CiCl) + i^[CI/\Cq - QACo]Q(r) - /i//(r)V (|CopV 
Here 



Q{r) = Jf{r')dr\ H{r) = J J ^J{r')dT\ 

0-^'^^* 



If K = kq = const then 

Q(r) = —(1 - e-^'^), H{t) = r— - \{l - e'^^^] 



One can see that Q{t) decays exponentially to a constant, whereas H{t) in- 
creases linearly with time, i.e. corresponding component is secular, and for 
stable solution it should become zero. Hence, we come to the following: taking 
into account ([HI), functions Co(r) and Wo(r) must obey the equations such 
that 



V(|CopV$o) = 0, (16) 
2 



A<Po = ^\Co\'. (17) 



Then equation f p^ reduces to 

A$i = ,2{CoCl + CiC*) + ^div[|Co|2veo] (1 - e-^"^) , (18) 

where Ci(r) can be obtained from stationary condition at the next order of 
expansion , and Oq = (i/2)ln(Co/Co) - is still an arbitrary function. The 
solution for Wi follows from the first equation in ( IT!?| l. 

The interpretation of obtained equations follows from the expression for the 
flow velocity in the first order 



Vi = vi + V2, vi = £V$o, V2 = 2£/(r)Veo(r) = 2£e-'^""Veo(r). 
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Right hand member in ( p^ is determined by the source of mass, which is 
associated with the second flow in the system. This flow tends to a fixed space 
distribution as r ^ oo. Equation ( 0^ is the law of conservation of mass for 
the flow Vi. Suppose sources of mass do not exist; then the second addend in 
right hand side of ( fTHl l becomes zero 

divflCopVGo] = 0. (19) 

If not, we must explicitly write down the source of mass by means of join of 
external and internal solution. The first flow is stationary and it is associated in 
first order with stationary fall in the field $o- Existence of dissipation leads to 
fall of particles with the fixed velocity vi instead of falling with the acceleration 
g = -V$o- 

5 Axial-symmetric solutions 

Our aim is to consider models with axial symmetry. Take cylindrical polar 
coordinates r, (p implying dependence functions of the system on r and z 

^^{r^ z) = u{r)h{z)^ "T^l^, z) = p{r)h{z). 

Thus one can find the following equations for w(r), p{r), h{z): 

u" lu' h" HP u' p' {h'Y fip ^ 
u r u h 2w' up 2u 

Separation of variables implies 

h{z) = hoe^^'. 

We should suppose dust density and potential to decay while moving off 
the disk in the line of z ^ +oo as well as z ^ — oo. Thus if A > 0, we have 
h{z) = /ige"'^'^' far from the disk surface z = 0. 
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So we obtain equations for u and p 



u" + -u' + \\ = ^p, (20) 
r 2 

u' p' ^0 lip , , 



up 2u 

Let us seek the solution for p as: p{r) = q{r)u'{r). Substituting p in above 
form into equation f i2T| l and using f i2n[ ) yield 



/ 1 2 

(? — q + i^q =0. 

r 

This equation is easy to solve and general solution is 
here Qo - integral constant. Finally we get 



where u{r) now follows the equation 



u" + , ^ u' + A^^ = 0. (23) 

r(r2 + Qo) 



Function p(r) follows, correspondingly, the equation 



P + / 2^n + A 1? = 0. 24 



At small Q or large r the equation for u{r) is similar to oscillation equation 
with wave number A and, hence, u{r) changes the sign quasi-periodically as 
well as its derivative. 

So we can establish behavior of dust distribution. The definition of p{r) 
implies p{r) > 0. The density becomes zero together with the gradient of 



potential according to f i22l ). Thus it follows that disk is partitioned on rings 
which are separated by thin gaps. Analysis of equations for flow shows that dust 
from certain ring does not penetrate the bounders. The gradients of potential 
in adjacent rings, however, have unlike signs due to its quasi-periodic behavior. 
But if we require continuity of density and its derivative at boundary points, 
we immediately obtain that density must attain a negative value along with 
u' . Therefore, the derivative of density has discontinuity at boundaries of rings 
and in each ring with constant sign of u' we should choose the sign and the 
magnitude of Q to satisfy requirement p{r) > 0. 

Generally boundary conditions for calculation of ring parameters come to 
continuity of potential and its derivative at ring boundaries (second derivative 
is discontinuous). It follows from equality of forces operating at disk boundaries. 

Let Tj, ? = 1, 2 ... be boundary points in which u'\r^ = 0. Then equations 
for our model are: 




Q. 



u.+\]uj = 0, r G h,rj+i] 



with boundary conditions 







Here aj is a constant surface density in j ring. 
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Fig.l. 1 - potential u{r) , 2 - density p{r). Ring parameters : A = 1, 
li = 2,Qo = 10, Qi = -19.85, Q2 = -19.7, Q3 = -155, (^4 = -108. 

Figure 1 illustrates the different possible solutions satisfying the boundary 
conditions. 



6 Discussion 

We have investigated the formation of rings with too narrow gaps between 
them as t — 00. And the width of gaps are considerably smaller then disk 
rings one. The masses of rings and distribution of density are determined by 
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the coefficient Q, that can be unique for each ring. Obtained solutions can be 
associated with dust distribution in real disk systems. However, to describe 
systems like internal Saturn's rings we must use other methods because our 
approach is not suitable for internal regions. Radial orbital flow induced by 
the central mass is the main component of dynamics for internal rings. Nearby 
the planet the radial velocity is great but we assume the mechanism of ring 
formation to be the same and it can be modified by corrections connected with 
the main orbital flow. And it will be the object of another paper. 
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